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Abstract 



With the early afterglow localizations of gamma-ray burst 
positions made by Swift, the clear delimitation of the prompt 
phase and the afterglow is not so obvious any more. It is 
important to see weather the two phases have the same origin 
or they stem from different parts of the progenitor system. We 
will combine the two kinds of gamma-ray burst data from the 
\ Swift-XRT instrument (windowed timing and photon counting 

modes) and from BAT. A thorough desription of the applied 
procedure is given. We apply various binning techniques to the 
different data: Bayes blocks, exponential binning and signal- 
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■ to-noise type of binning. We present a handful of flux curves 

and some possible applications. 



1 Introduction 

Gamma-ray bursts are the most energetic phenomenon in the Uni- 
verse. After their discovery (Klebesadel et al., 1973), and the detec- 
tion of the first counterpart in other wavelengths than gamma now 
the Swift satellite (Gehrels et al. 2004) is almost routinely observing 
X-ray afterglows. In the study of GRB prompt emission and after- 



glow, it is a straightforward idea to combine flux curves of GRBs from 
gamma-ray and X-ray data. A wide spectral coverage leads to a more 
complete picture of the phenomenon. In our analysis have extrapo- 
lated the gamma-ray flux into the X-ray band to have a commond 
ground for analysis. It can be done the other way around, extrapo- 
lating the X-ray flux to the BAT energy range. Most bursts in our 
sample come from the long and possibly the intermediate duration 
group (Horvath 1998, Balazs et al. 1998, Balazs et al. 1999, Horvath 
2002, Horvath et al. 2008, Ripa et al. 2009) 



2 Data reduction and binning 

Swift gathers 7-ray and X-ray data relevant for our analysis. We 
choose three different approaches to bin the flux curves: Bayesian 
method for the gamma-ray data, equal binning in logarithmic coor- 
dinates in the case of the windowed timing (WT) XRT data and a 
signal-to-noise type of binning in the case of photon counting (PC) 
XRT data. The flux curves and the spectra were generated using stan- 
dard HEASof t tools and the most recent calibration database. Initial 
calibration was made using xrtpipeline and batgrbproduct pipeline 
scripts with the latest calibration databases. 

In gamma-rays the results of the primary pipeline processing con- 
tain, among others, a 64 ms resolution gamma-ray flux curve. This 
was used as the input to the Bayesian block analysis. We cut our 
combined 15 — 150 keV flux curve using Bayes blocks as presented in 
Scargle (1998). We set a large prior for the algorithm so it will stop at 
an early point (which corresponds to a small number of change points) 
making sure that we have enough resolution for each of the bins. After 
deducing the time intervals we have used batbinevt to bin the data 
and get a raw spectrum (pha) file. The further steps recommended 
in the BAT analysis thread were also carried out. For each interval, 
we created the appropriate response matrices and fitted both a power 
law and a power law with a high-energy cutoff. We use the criterion 
from Sakamoto et al. (2008) to choose between the two models. If 
the x 2 improves by more than 6 by using a cutoff power law, we use 
the latter instead of the simple power law model. The next step is to 
extrapolate the model from the gamma-ray band (15 — 150 keV) into 
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the X-ray band (0.3 - 10 keV). 

In X-rays the WT mode is active when the count rate of the source 
is high (over ~ 10 counts/s). This means we have a good signal-to- 
noise ratio and we can bin the counts in equal bins in logarithmic 
space. We have fitted a spectrum for the whole duration of the WT 
mode and got a conversion factor from rates to flux. There is a more 
detailed description of the procedure in the next section. 

For the PC mode we bin our data to have a signal-to-noise ratio of 
at least 3. We do this by incrementing the endpoint of our interval in 
time until the count rate reaches the required level. At this point we 
store this interval and repeat the procedure until the end of the ob- 
serving period. We correct for the pile-up in the detector as described 
in Vaughan et al. (2006) for the PC mode. For the WT mode pileup 
correction is made according to Romano et al. (2008) 

The next step is to convert the count rates to flux. To do this, we 
divide the cumulative flux curve in n parts, each with equal number 
of counts, n is chosen by hand depending on the intensity of the 
afterglow from 2 to 6. For each time slice we fit a spectrum to get the 
count equivalent in erg/cm 2 . The spectra are fitted with Xspec using 
individual anciliary response functions and the most recent response 
function available. We used an absorbed power law model of the form: 

(wabsc) x (zwabss) x (powerlaw). 

The absorbing column density (NH) of the Milky Way was taken from 
Dickey et Lockman (1990) (denoted here by wabs G ), and, where red- 
shift was known, the source absorbtion was also fitted (denoted here 
by zwabsg). If the redshift was unknown zwabs^ was substituted with 
a simple absorbing wabss component. The spectra were binned using 
grppha so all channels had minimum 20 counts. 

Evans et al. (2007) use the same method, but they get the con- 
version factor (erg/cm 2 /s equivalent of 1 count/s ) by integrating over 
the entire spectrum (corresponding to n — 1). We report here that 
our conversion factors are in good agreement with those in Evans et 
al. (2007) and their related web-page. 
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Figure 1 : Photon indices from the late prompt emission phase and the 
early afterglow phase. The continuous line represents equality between 
the two parameters 

3 Interesting cases 

Here we present several GRBs to illustrate our method of binning 
and combining data. For now, we include only bursts from the long 
and intermediate duration population. The behaviour of the afterglow 
requires the use of a log scale but several GRBs have significant flux 
before the trigger (precursors with t < T trigger = 0). For this reason 
we choose to include a shift in the time axis where necessary. 

Another way of examining the problem is to consider gamma-ray 
and X-ray photon indices. If the emission stems from the same pop- 
ulation of electrons, we would expect the two types of indices to be 
equal within errors. The measurements should be carried out roughly 
at same time. We have taken a sample 10 bright Swift bursts and com- 
pared the two indices. We found that there was a visible discrepancy. 
(Figure □]) 



4 



4 Discussion and Conclusion 



We have selected a few bright GRBs (Figure H]) to present our binning 
method for data from two bands. There are cases where the possible 
extrapolation of the WT mode data seems to match the gamma-ray 
flux curve (for the first three GRBs: 060210, 060418 and 060502). 
This could be proof of a connection between the processes which gov- 
ern the prompt phase and the processes of the early afterglow. It 
looks possible that the steep decline is indeed a sequel of the prompt 
phase seen in X-rays. No such claim can be made for the following two 
events (060604 and 060908) because of the scarceness of the data. At 
the last two events (061121 and 060124) however, there is an appar- 
ent discrepancy between the gamma-ray curve and the X-ray curve. 
This discrepancy is best seen at GRB061121. This could mean the 
extrapolation of the gamma-ray spectrum was not adequate. 
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GRB060124 - Combined flux curve 



GRB060210 - Combined flux curve 
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Figure 2: A few examples of combined flux curves of Swift GRBs. 
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